Flat Central Density Profiles from Scalar Field Dark Matter Halos 
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Scalar fields endowed with a cosh potential behaves in the linear regime, exactly as the cold dark 
matter (CDM) model. Thus, the scalar field dark matter (SFDM) hypothesis predicts the same 
structure formation as the CDM model. The free parameters of the SFDM model are determined 
by cosmological observations. In a previous work we showed that if we use such parameters, the 
scalar field collapses forming stable objects with a mass around 1O 12 M0. In the present work we 
use analytical solutions of the flat and weak field limit of the Einstein- Klein- Gordon equations and 
show that the SFDM density profile corresponds to a halo with an almost flat central density and 
that it coincides with the CDM model in a broad outer region. This result could solve the problem 
of the density cusp DM halo in galaxies without any additional hypothesis, supporting the viability 
of the SFDM model. 

PACS numbers: 



I. INTRODUCTION 

The Lambda Cold Dark Matter (ACDM) model has 
recently shown an enormous predictive power. It can ex- 
plain the structure formation of the Universe, its acceler- 
ated expansion, the micro Kelvin fluctuation of the Cos- 
mic Microwave Background Radiation, etc. Nevertheless, 
some issues around this model related to the formation of 
galaxies have arisen since the time it was originally pro- 
posed and remain to date. The CDM paradigm predicts 
a density profile which corresponds to the Navarro- Frenk- 
White (NFW) profile, Navarro et al, 1997, given by 
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However, this profile seems to have some differences to 
the observed profiles of LSB galaxies. In this work we 
show that a flat central profile naturally arises within 
the scalar field dark matter hypothesis, implying that the 
central region of galaxies can distinguish between CDM 
and SFDM. 

We work within the specific context of the so-called 
'strong, self-interacting scalar field dark matter' (SFDM) 
hypothesis that has been developed by several authors, 
Guzman and Matos 2000, Matos and Urena-Lopez 2000, 
2001; Urena-Lopez, Matos and Becerril 2002; Matos and 
Guzman 2001; Alcubierre et al. 2002, 2003; Urena-Lopez 
2002; Bohmer and Harko 2007; (see also Peebles 2000). 
A first proposal of the SFDM hypothesis appeared in a 
couple of papers by Ji and Sin 1994. They took a mas- 
sive scalar field and were able to fit observations comming 
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from some galaxies, taken into account also the contri- 
bution of baryons. A key point in this work was the use 
of the so-called excited configurations in which the radial 
profile of the scalar fields has nodes. From that, they 
determined that the mass of the scalar field should be of 
order of ~ 10~ 24 eV. 

A next proposal appeared in a paper by Schunck 1998. 
He shows that a massless complex scalar field can be used 
as a dark matter model in galaxies to fit the rotation 
curves. In this model the internal frequency of the field 
plays the role of an adjustable parameter, and the radial 
profile of the scalar field also has nodes. 

However, as pointed out in Guzman and Urena-Lopez 
2005 (see also Guzman and Urena-Lopez 2004), this last 
proposal cannot be realistic because a massless scalar 
field (whether real or complex as in Schunck 1998) can- 
not form a gravitationally bound configuration (see Seidel 
and Suen 1994). 

The key idea of the SFDM scenario is that the dark 
matter responsible for structure formation in the Uni- 
verse is a real scalar field, $, minimally coupled to Ein- 
stein gravity with self-interaction parametrized by a po- 
tential energy of the form (see also Sahni and Wang 2000) 



V($) = V [cosh(Av^u"$) - 1] . 



(2) 



where Vq and A are the only two free parameters of 
the model, kq = 8ttG and we employ natural units 
h = c = 1. The effective mass of the scalar field is given 
by m| = KoVbA 2 . 

The advantage of the SFDM model is that it is insen- 
sitive to initial conditions and the scalar field behaves as 
CDM once it begins to oscillate around the minimum of 
its potential. In this case, it can be shown (see Matos and 
Urena-Lopez 2000, 2001) that the SFDM model is able to 
reproduce all the successes of the standard ACDM model 
above galactic scales. 
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Furthermore, it predicts a sharp cut-off in the mass 
power spectrum due to its quadratic nature, thus ex- 
plaining the observed dearth of dwarf galaxies, in con- 
trast with the possible excess predicted by high resolu- 
tion N-body simulations with standard CDM, see Matos 
and Ureha-Lopez 2001. 

The best-fit model to the cosmological data can be de- 
duced from the current densities of dark matter and ra- 
diation in the Universe and from the cut-off in the mass 
power spectrum that constrains the number of dwarf 
galaxies in clusters. The favored values for the two free 
parameters of the scalar field potential are found to be, 
Matos and Ureha-Lopez 2001: 

A ~ 20, 

V a ~ (3xlO- 27 m Pi ) 4 , (3) 

where mpi = G -1 / 2 » 10 _5 g is the Planck mass. 

This implies that the effective mass of the scalar held 
should be m$ ~ 9.1 x 10~ 52 m P i = 1.1 x 10~ 23 eV. 

Let us explain why we suspect that the scalar held 
could be the dark matter at galactic scales as well. There 
are three main reasons. 

The hrst reason is that numerical simulations suggest 
that the critical mass for the case considered here, us- 
ing the scalar potential (|2|) , and the parameters given by 
Eq. |3j), is approximately, Alcubierre et. al. 2002 

o 

Mcnt ~ 0.1^4= = 2.5 x 10 13 M o . (4) 

This was a surprising result. The critical mass of the 
model shown in Matos and Ureha-Lopez 2000, 2001, is of 
the same order of magnitude of the dark matter content 
of a standard galactic halo. Observe that the parame- 
ters of the model, Eq. (J3J) , were fixed using cosmological 
observations. The surprising result consisted in the fact 
that using the same scalar field for explaining the dark 
matter at cosmological scales, it will always collapse with 
a preferred mass which corresponds to the halo of a real 
galaxy. Thus, this result is a prediction of the cosmolog- 
ical SFDM model for galaxy formation. 

The second reason is that during the linear regime 
of cosmological fluctuations, the scalar field and a dust 
fluid, like CDM, behave in the same way. The density 
contrast in CDM and in the SFDM models evolve in ex- 
actly the same form and then both models predict the 
same large scale structure formation in the Universe (see 
Matos and Ureha-Lopez 2001). The differences between 
the CDM and SFDM models begin to appear in the non 
linear regime of structure formation, so that there will be 
differences in their predictions on galaxy formation. 

The third reason is the topic of this work. A scalar 
field object (e.g. an oscillaton) contains a flat central 
density profile, as seems to be the case in galaxies. 

In the case of the SFDM, the strong self-interaction 
of the scalar field results in the formation of solitonic 
objects called 'oscillatons', which have a mass of the order 



of a galaxy (see for example Urena-Lopez 2002, Ureha- 
Lopez et al 2002 and Alcubierre et al 2003. Also Seidel 
and Suen 1991, 1994, Hawley and Choptuik 2000, and 
Honda and Choptuik 2001). In this work we will show 
that these models contain an almost flat central density 
profile, i.e., they do not exhibit the cusp density profiles 
characteristic of the standard CDM hypothesis. 

Before starting with the description, we want to emp- 
hazise the fact that the scalar field has no interaction with 
the rest of the matter, thus, it does not follow the stan- 
dard lines of reasoning for the particle-like candidates for 
dark matter. The scalar field was not thermalized, that 
is, the scalar field forms a Bose condensate, and thus 
behaves strictly as cold dark matter from the beginning. 

The rest of the paper is organized as follows. In the 
next section we use the fact that Galaxies have a weak 
gravitational field and thus their space-time is almost 
flat. The main goal of this section is to study the physics 
provoking the flatness behavior of the density profiles at 
the center of the oscillatons. Some results of this sec- 
tion intersect with those presented by Jin and Sin 1994, 
where they studied the behavior of the weak field limit 
of a complex scalar field. We remark that we do our 
analysis for a real scalar field and some differences do 
arise due to a different current conservation. In section 
IIIII we present our contribution making an analysis of 
the Einstein-Klein-Gordon (EKG) equations in the rela- 
tivistic weak field limit and solve them for the perturbed 
metric coefficients. Then, we compare these solutions 
with the ones obtained by solving numerically the com- 
plete EKG system and using the whole potential |(3J), wc 
show that the relativistic weak field limit is indeed a very 
good approximation. In section IIVI we compute the en- 
ergy density of the scalar held obtained in section uTTl and 
compare it with actual observations of LSB galaxies from 
which the density is inferred from the rotational curves, 
showing a good matching in the external regions and a 
match at least similar in some of the internal regions, in 
any case, better than a fit with a cusp-like behavior of 
the density. Finally in section [V] we give our conclusions. 



II. PHYSICS OF THE SCALAR FIELD. FLAT 
SPACE-TIME CASE 

In this section we derive the physics of the scalar field 
in an analogous way as it was done for the complex scalar 
field by Jin and Sin 1994 (see also Lee and Koh 1996). 
In a normal dust collapse, as for example in CDM, there 
is in principle nothing to avoid that the dust matter col- 
lapses all the time. There is only a radial gravitational 
force that provokes the collapse, and to stop it, one needs 
to invoke some virialization phenomenon. In the scalar 
field paradigm the collapse is different. The energy mo- 
mentum tensor of the scalar held is 

= - 9 -f- + 2U($)] . (5) 
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We will consider spherical symmetry, and work with the 
line element 



ds 2 = 



^dr 2 



r 2 dn 2 



(6) 



with \i = n(r,t) and v = v{r,t), being this last function 
the Newtonian potential. The energy momentum tensor 
of the scalar field has then the components 
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(10) 

and also T 3 3 = T 2 2- These different components are 
identified as the energy density the momentum den- 
sity V$>i the radial pressure p r and the angular pressure 
p± . The integrated mass is defined by 



M(x) 



4tt f p< s> (X)X 2 dX. 
Jo 



(11) 



The radial and angular pressures are two natural com- 
ponents of the scalar field which stop the collapse, avoid- 
ing the cusp density profiles in the centers of the collapsed 
objects. This is the main difference between the normal 
dust collapse and the SFDM one. The pressures play 
an important roll in the SFDM equilibrium. In order to 
see this, and considering that galaxies are almost flat, 
we conclude that the Newtonian approximation should 
be sufficient to describe the processes. In this section we 
will take the flat space-time approximation. 

Thus, we study a massive oscillaton without self- 
interaction (i.e. with potential V = |m|$ 2 ), in the 
Minkowski background (/x ~ v ~ 0). Even though 
it is not a solution to the Einstein equations as we 
are neglecting the gravitational force provoked by the 
scalar field, the solution is analytic and it helps us to 
understand some features that appear in the non-flat 
oscillatons. 

In a spherically symmetric space-time, the Klein- 
Gordon equation r) a °d a d/3& — dV/d§ = 0, (where 
rj a Pd a dp stands for the D'Alambertian), reads 



(12) 



where over-dot denotes d/dt and prime denotes d/dr. 
The exact general solution for the scalar field $ is 



*(*,r) 



(13) 



and we obtain the dispersion relation k 2 = ui 2 — m\. 
For ui > m$ the solution is non-singular and vanishes at 



infinity. We will restrict ourselves to this case. It is more 
convenient to use trigonometric functions and to write 
the particular solution in the form 



$(t, x) = <&o — — cos(wi), 



(14) 



where x = kr. It oscillates in harmonic manner in time. 
The scalar field can be considered to be confined to a fi- 
nite region, see Urena-Lopez 2002, and Urena-Lopez et al 
2002. The analytic expression for the scalar field energy 
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FIG. 1: The energy density profile for the Scalar Field Dark 
Matter model (see also Schunck 1997). The parameters used 
in this plot are m* = 20, u = 20.001, $g = 5 10" 2 . The 
density is given in arbitrary units and r is given in kpc. 



density derived from Eq. (fl"4|) is 



<fo 2 k 2 
2x 2 



'sin(x) 



cos(a;) 



(15) 



-k 2 sin 2 (a;)) cos 2 (u>t) + u> 2 sin 2 (x)] 



which oscillates with a frequency 2u>t. Observe that close 
to the central regions of the object, the density of the 
oscillaton behaves like 



1 



p<s> 



$ 2 fc 2 [w 2 - k 2 cos 2 (cj t)] + 0(x 2 ), 



(16) 



which implies that when x — > the central density oscil- 
lates around a fixed value. 

On the other hand, the asymptotic behavior when 
x — * oo, is such that /?$ ~ 1/x 2 , i.e. far away from the 
center, in this approximation, the flat oscillaton density 
profile behaves like the isothermal one. The mass func- 
tion oscillates around M usual for the galactic 
halos. 

In order to understand what is happening within the 
object, observe that the KG equation can be rewritten in 
a more convenient form in terms of the energy density, 
as 



dt 



r 2 dr 



2 P$) = 0. 



(17) 



4 



This last equation has a clear interpretation: Since its 
form looks like the conservation equation, p + V • J = 0, 
equation (fT7|) represents the conservation of the scalar 
field energy. It also tells us that there is a scalar field 
current given by 



J<$> = —Vq>r 

, 9 few r , . . , > , sin(x) sin(2u;£) _ 
= $0— [xcos(x) - sm(x)\ -r. 



Observe that the quantity involved in this current is 
the scalar field momentum density ([8]). Although the 
flux of scalar radiation at large distances does not van- 
ish, there is not a net flux of energy, as it can be seen by 
averaging the scalar current on a period of a scalar oscil- 
lation. We also see that the only transformation process 
is that of the scalar field energy density into the momen- 
tum density, and viceversa. For the realistic values ([3]) 
this transfer is very small. 

In Fig. [1] we show the behavior of the SFDM den- 
sity profile for a typical galaxy and in Fig. [2] we show 
the comparison between the NFW, the isothermal and 
the SFDM density profiles for the same galaxy. Observe 
that the SFDM and NFW profiles remain very similar 
up to 10 kpc, then the SFDM profile starts to follow the 
isothermal one. 



III. WEAK FIELD LIMIT EQUATIONS 

In this section we derive a novel method for integrat- 
ing the perturbed EKG equations and show that the so- 
lutions are in very good agreement with the numerical 
ones. This allow us to use these solutions to fit the rota- 
tion curves of several observed LSB galaxies. 

Within general relativity, the evolution of the scalar 
and gravitational fields are governed by the coupled EKG 
equations, the last one appearing from the conservation 
of the energy-momentum tensor 



Raf3 — KoS a f3, 

T af3 - P = $-"(V 2 



TO 2 )$ = 0, 



(18) 
(19) 



- (l/2)g aP T\, 
the covariant 



here R a p is the Ricci tensor, S a p = T a p 

and V 2 = (1/V=9)^[V=9<7 M! U] is 
D'Alambertian operator. 

For simplicity, we continue to consider the non-static 
spherically symmetric case, given by Eq. ©. As usual in 
the weak field limit, we suppose the metric to be close to 
the Minkowski metric r\ a p 



9ap = rj aj 3 + h a j3, 



(20) 



where \h a p\ <C 1, then we will consider an expansion of 
the functions in the metric of the form 



S, -4 



Isothermal Profile 
NFW Profile 
SFDM Profile 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 

log(r) (kpc) 

FIG. 2: Comparison between the energy density profile for 
the Scalar Field Dark Matter model with the NFW and the 
Isothermal models. The parameters for the isothermal model 
are pi so = 0.3/ (r 2 +8 2 ) and for the NFW profile are Pnfw = 
10/ (r(r + 8) 2 ). The parameter used for the SFDM model are 
the same as in the previous figure. 



In the next section we will see that if the gravitational 
force is taken into account, the oscillaton is more 
confined (see also Urena-Lopez 2002, Urena-Lopez et al 
2002, and Alcubierre et al 2003). The parameters used 
in the figures, correspond to a middle size galaxy. 



M*,r) = l + e 2 2V{t,r) + 0(e 4 ), 
,2M*,r) = l + e 2 2U(t,r) + 0(e 4 ) 1 



(21) 



where e is an expansion parameter. We also consider 
that the spatial and time derivatives of the geometric 
quantities are regarded like 



d_ d_ 

dt dr 



(22) 



then to first order in e 2 the Ricci tensor components are 
respectively 

Rtt = [V, rr -U >tt +-V, r }, 



Rrr = [U,tt~V, rr +-U, r ], 

r 

Rtr = [~U, t ]. 
r 



(23) 



On the other hand, the source is computed in the flat 
space in this case as well, thus the scalar field satisfies to 
e 2 order the Eq. (fT2f as in the previous section. 

It is important to emphasize that the relation |22|) is 
the lowest one in the geometric fields, but it does not 
consider small velocities for the sources. This is different 
from the Newtonian limit where the derivative relation 
for the scalar field is d r ~ ed t and d t ~ ed r for the 
geometric fields (see Seidel & Suen 1990, and Guzman & 
Urena-Lopez 2004). 
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Consistent with the T a p computed in the flat space- 
time, the right hand side elements in Einstein's equations, 
are written as 

S tt = $, t $, t -im 2 |$| 2 , 
1 



A. Scaling properties 



-m 2 |$| 2 , 
2 1 1 



S rr = 

S tr = (24) 

In this case it is convenient to introduce the dimensionless 
quantities 

x = mr, t = mt, fi = — , (25) 

m 

where we note that the bosonic mass m is the natural 
scale for time and distance. In terms of these new vari- 
ables the general solution to Eq. (fT5)) takes the form 



$(t, x) = — exp(±ixy/ ft 2 — 1) exp(±iOr). 



(26) 



The physical properties of the solution depend on the 
ratio fl = ui/m. For < 1 the solution decays expo- 
nentially but it is singular at r = 0. On the other hand, 
fl > 1 allows for non-singular solutions which vanish at 
infinity. We will restrict ourselves to this case. We will 
write the particular solution in the form 

(27) 



(28) 



y / Ko<f>(r, x) = fa(x) cos(Or), 
where the spatial function is given by 



(i) = fa 



suv(x^n 2 - 1) 



Because of the functional form of the scalar field ([27] 
we introduce the following ansatz for the metric pertur- 
bations |[2"T]) 



V(t,x) = V (x) + V r 2 (x)cos(2rJr), 
U{t,x) = U {x) + U 2 (x)cos(2CIt). 



(29) 



We adopt the following method to solve the EKG equa- 
tions. In the first approximation we substitute the solu- 
tion of a lower approximation (the flat case) into the next 
approximation (with e 2 ) and solve the resulting differen- 
tial equations. As we will see, this standard approxima- 
tion works well in our case. In terms of these expressions 
the Einstein's equations R a p = K[)S a p finally read 



V , 



4 

~U 2 

x 



-V , x 



(30) 



v 2 , xx +-v 2 , x +iwu 2 

X 



X 



-v 2 , xx +-u 2 , x -awu 2 

x 



l 

2 
1 

+ 2 
-V 

-V 



From system (|30| we know that the scalar field's maxi- 
mum amplitude fa(0) = faoV^ 2 — 1 could be taken as the 
expansion parameter e and in this case fl must be of or- 
der 1. Then it is always possible to solve the system (|3"0")) 
ignoring e and replacing fa by its normalized function 



fa(x) 



sin(:r\/fF— T) 
XyJQ 2 - 1~ ' 



(31) 



Solutions fa, Uq, U 2 , Vq, V 2 of this normalized system 
depend only of the arbitrary characteristic frequency Q 
which modulate the wave length of fa. On the other hand, 
for each value of f2 there is a complete family of solutions 
of the scalar field fa and the metric perturbations h a p 
which are related to each other by a scaling transforma- 
tion characterized by fao 



y/K0$ = 

htt = 



fao\/n 2 — l^cos(fir), 

fal{n 2 - 1)[2U + 2C7 2 cos(2Qt)], 

-fal(n 2 - 1)[2V Q + 2F 2 cos(20r)], 



(32) 



In this context the weak field limit condition h a p <C 1 
translates into 



fa 2 a {n 2 - 1)\2V\ < 1, & fal{n 2 - l)\2U\ < 1. (33) 



Here we will introduce a specific notation for the spatial 
functions of the metric perturbations: 





fa 2 (n 2 - 


l)2U , 




fa 2 (n 2 - 


l)2C/ 2 , 


<ht — 


fa 2 (n 2 - 


1)2V , 


'Hi — 


fa 2 (n 2 - 


1)2V 2 . 



(34) 



B. Metric perturbations solutions 



The system of equations Eqs. (I30|) can be solved and 
the spatial functions of the metric perturbations have 
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analytic solutions given by 



U 2 = 



8(n 2 - 1) 



sm 2 (x\/tt 2 - 1) 



V^ 2 - 1 sin(22V^ 2 - 1) 



Vn = 



(2fl 2 - 1) rsin(2x\/f7 2 - 1) 
8(ft 2 - 1) I 2xVW~ r T 
Ci{2x^/n 2 - 1) + 



Cvoi 



CV02, 



V 2 = 



U 



ln(2x\/n 2 - 1) 

J J/ 

1 rV^ 2 - 1 sin(2a;V^ 2 - 1) 
8(Q 2 - 1) i 2 x 

Ci(2xy/n 2 - 1) - 
ln(2a;v / ^ 
1 



+ CV22, 
1 1 cos(2x%/^ 2 - 1) 



8(Q 2 - 1 

1 sin(22Vft 2 - 1) 
2 



2x 2 2 x : 
CV01 



1 



+ CW, (35) 



Unique solutions for the EKG system will be obtained 
fixing the </>o and Q parameters and the constants CV02 
and Cyoi within the validity range of the approxima- 
tion. As it is known the potentials measurement does not 
have physical sense by themselves, then unique solutions 
will be determined through metric dependent observable 
quantities. Using the expressions given by Eqs. (|34[) . we 
can obtain the perturbed metric functions in terms of 
these solutions. In Fig. [3] we present a plot of these met- 
ric perturbation functions, as well as of the scalar field, 
for two values of SI. It is important to mention that 



£2=1.4 




where Ci is the cosine integral function and Cvoi, CV02, 
CV22, and Cum are integration constants. 



C. Weak Field Validity Range 

From equation (|13| it is evident that, in the limit in 
which we are working with, the KG equation is decoupled 
from Einstein equations. Imposing regularity at the ori- 
gin and asymptotic flatness to the KG solution we have 
chosen (|2"T| with (|2"8")l as our scalar field particular solu- 
tion where <po an d SI > 1 are still free parameters. 

On the other hand, regularity at x = requires 
h rr (x = 0, r) = which implies 

si 2 

Cvoi = 0, C U01 = 2 - — , (36) 

then for the perturbations CV02 and CV22 are still free 
integration constants. Now we will describe the asymp- 
totic behavior of these perturbations. Due to U 2 being 
at least one order of magnitude smaller than U$ and its 
value oscillating around zero, it is Uo which determines 
the behavior of h rr . The Uo value starts to oscillate, 
very near to the origin, around Cum keeping this behav- 
ior asymptotically. Then the asymptotic value of h rr is 
the finite Cum value. Contrary to this htt, due to the 
logarithm terms in Vq and V2, h u , is singular at infin- 
ity. Thus, the weak field condition (|33[) is fulfilled only 
in a finite spatial region around the origin, i.e., due the 
the approximation the solution is contained in a box, for 
which the walls are sufficiently far away from the center 
of the solution. We will say that this is the region where 
our weak field approximation is valid. 
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FIG. 3: Profiles of <f>{x), h?){x), hf) {x),h { $ (x) and h$ {x) 
with (j>o = 1; see text for details. 

the width of the validity region, where (|3"3"|) is fulfilled, 
depends on Q and <fio. This is because it is the factor 
0o which modulate the perturbations, (see (|32|) ). What 
it is not evident until the solutions ([3"5]) are observed is 
that the £1 value, independently from (f>o, could make the 
validity range width bigger. This is beca use in t he loga- 
rithm argument there is the expression \/Vl 2 — 1, then as 
O is closer to one the logarithm terms rise more slowly. 

The order of magnitude for the other parameter (j) 
in the metric perturbations, can be naturally determined 



7 




£1=1.181 Analytical 

+ Numerical 
— - Q=1.434 Analytical 

X Numerical 



"h (t , £1=1.181 
-h it 12 '. £1=1.181 

Numerical 
V, 0=1.434 
- h B ,Q=1.434 

Numerical 



plots. As already was noticed, the value of character- 
izes each family of solutions. Mainly determines the 



wave length of <f> and the increase rate of h^' and h)£ ' ; as 



is closer to 1 this rate is smaller, 
are shown in Fig. [3j 



These characteristics 



p, t„=5.222e-3 

+ Numerical 

p*5, ^ -1.591e-3 
x Numerical 




> p-p N . «> =5.222e-3 

< (P-P„)»100. f, () =1.591e-3 



FIG. 4: (Left) Spatial function </> of the scalar field <F 
Solid a nd dash lines are analytical solutions with <j)p — 
0.001/-y/(fi 2 - 1) = 0.001591, and <f) = 0.002/,/(fi 2 - 1) = 
0.001946, respectively. The crosses are the corresponding 
numerical solutions with 0i (0) = 0.001, Q = 1.181008 and 
4>i (0) = 0.002, fi = 1.433822; see text for details. (Right) 



Spatial functions h$ and h\l' of the metric perturbations. 



,( 2 > 



from the asymptotic value taken by h rr , which is reached 
very close to the origin 



hm \h rr \ = -^n 2 . 



(37) 



As 51 is nearly 1 the h rr magnitude is given by (f>Q. It 
is well known that for week field systems like our Solar 
System the metric perturbations go like h a p ~ 10~ 6 . 
This value restricts our <p to be <po < 10~ 3 . 



D. Analytical Solutions vs Numerical Solutions 



Analytic solutions </>o, h% (x), hrr(x), h^' {x) and 
h$ (x) are shown in Fig. [3] The value of 4>q is 1 in both 



(0), 



,(2) 



FIG. 5: (Left) Energy-Momentum density profiles for two 
scalar field configurations one with <j>o — 1.1591 10~ 3 and 
Q = 1.181008 and the other with O = 5.2225 10~ 3 and fi = 
1.153246. (Right) Difference between density profiles from 
the exact EKG equations (pjv) and the density computed in 
the flat space-time (p); see text for details. 



The exact EKG equations in spherical symmetry with 
a quadratic potential, were solved numerically in Ureha- 
Lopez 2002, and Ureha-Lopez et al 2002) and found the 
so called oscillatons. In those works boundary conditions 
are determined by requiring non-singular and asymptot- 
ically flat solutions, for which the EKG become an eigen- 
value problem. The free eigenvalue is the scalar field's 
central value (f>i(x = 0) which labels the particular equi- 
librium configuration, and the fundamental frequency 
f2 is an output value. In those works it is also noted 
that weak gravity fields are produced by oscillatons with 
(j)i(x = 0) <C 1. In Fig. d] we compare some of these nu- 
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merical solutions (NS) with the analytical solutions (AS) 
within a central region. The constant values of the AS 
are fixed to fit better the NS inside the weak field validity 
range. From these plots we can conclude that our solu- 
tions are a very good approximation for the exact EKG 
equations in the weak field limit. The principal advan- 
tage of this approximation is the analytical description 
of the solutions. 



IV. SCALAR FIELD AS DARK MATTER: HALO 
DENSITY PROFILE 

In this section we explore whether or not the scalar 
field could account as the galactic DM halos. Specifi- 
cally we compare the SFDM model density profile and 
the profiles inferred throughout the rotation curves pro- 
files of galaxies which are mostly formed by DM. 

As long as we are concerned with perturbations of the 
fiat space-time due to the scalar field, we do not consider 
the baryonic matter gravitational effects, we expect that 
our approximation will be better suited for galaxies with 
very small baryonic component. 

We will compare the energy-momentum density for the 
scalar field given by Eq. ([7]), in the relativistic weak 
limit approximation, where for the metric functions, 
Eg . (|2"Tj) . we use the solutions to the perturbations given 
by Eqs. (|29l35p . This is consistent with the fact that 



FIG. 6: Density profile fits for the galaxies ESO0140040, 
ESO2060140, the right panel is a zoom of its central region, 
it is plotted in order to visualize complete error bars. The 
horizontal line is in kpc, vertical line in M@/pc 3 . See text 
and Table [J for fit details. 



the difference between the density from the complete 
EKG equations and from our approximation is smaller. 



ESO3020120 ESO3020120 




FIG. 7: The same as Fig. \§\ but for the galaxies ESO3020120, 
ESO3050090. 

The density profile fits allow us to obtain an estimation 
of the parameters at the galactic level: the fundamental 
frequency f2 and the scalar field constant <pQ. The third 
parameter involved in the density profiles is the scalar 
field mass, we will fix it to be m = 10~ 23 eV. This value 
was fitted for the SFDM model from cosmological obser- 
vations in Matos & Ureha-Lopez 2001. 




FIG. 8: The same as Fig. [6] but for the galaxies ESO4250180, 
ESO1200211. 



gravity does not modify the scalar field behavior. This 
approximation in the weak gravitational field limit is very 
good as we can see in Fig. [5] In those plots we show 
the energy-momentum density profiles for two scalar field 
configurations with different maximum amplitudes at the 
origin 0(0) = 4>q\JVL 2 — 1. It is important to notice that 
as 4>(0) decreases the gravitational field gets weaker, and 



A. Density Profile fits 

The first qualitative feature of the energy-density pro- 
file that we want to emphasize is that in the central region 
it is non cuspy (see Fig. [5|). It is important to take into 
account that instead of density profiles, rotation curves 
are the direct observable for galaxies. Nevertheless, for 
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galaxies dominated by DM, their rotation curves could 
model the DM density profile more trustfully. We choose 
a subset of galaxies from the set presented in McGaugh 
2001, the common characteristic for the selected galaxies 
is that the luminous matter velocity contribution to the 
rotation curves is almost null. 

With the scalar field mass m fixed, the profiles fits were 
made for the f2 and (f>o values with good x 2 statistic, see 
Table HI In most of the cases the non central observa- 
tional data were the good fitted points, those data points 
also have smaller error bars. The density profile fits are 
in Figures I6I7I8I91 where we show several galaxies with 
the density computed from the observed rotational profile 
versus the density obtained with our SFDM description. 
In some of them we were able to compare regions within 
less than 0.5 Kpc. 

In Table U the fundamental frequency is listed for 
each galaxy. We found that the temporal dependence for 
the energy-momentum density profile is harmonic with 
a temporal period T = 7r/f2. The column Ap(0) cor- 
responds to the maximum change in the central density 
for a period of time T. Finally for all the galaxies the 
(j>a value is well inside of the weak gravitation field limit 
0o < 1(T 3 . 

V. CONCLUSIONS AND FUTURE PROSPECTS 

We have found analytic solutions for the EKG equa- 
tions, for the case when the scalar field is consider as 
a test field in a Minkowski background, and in the rel- 
ativistic weak gravitational field limit at first order in 
the metric perturbations. With these solutions we have 
shown that non-trivial local behavior of the scalar field 
holds the collapse of an object formed from scalar field 
matter. The scalar field contains non trivial, natural ef- 
fective pressures which stop the collapse and prevent the 
centers of these objects from having cusp-Ike density pro- 
files. Even within this simple approximation it has been 
possible to fit, with relative success, the density profiles 
for some galaxies showing non cuspy profiles. 

Together, all the features of the SFDM model allow one 
to consider this model as a robust alternative candidate 
to be the dark matter of the Universe, as was suggested 
by Guzman and Matos 2000, Matos and Guzman 2000, 
2001, and Matos et al 2000. Furthermore, it has been 
shown previously that dark halos of galaxies could be 
scalar solitonic objects, even in the presence of baryonic 
matter, Hu et al 2000; Lee and Koh 1996; Arbey et al 
2001, 2002; Sin 1994; and Ji and Sin 1994. Actually, the 
boson mass estimated in all these different approaches 
roughly coincides with the value m$ ~ 10~ 23 eV, even 
if the later was estimated from a cosmological point of 
view, Matos and Urena-Lopez 2001. We can appreciate 
the non-trivial characteristics of the proposed potential 
@: Its strong self-interaction provides a reliable cosmo- 
logical scenario, while at the same time it has the desired 
properties of a quadratic potential. Finally, the results 



presented here fill the gap between the successes at cos- 
mological and galactic levels. 



ESO 18705 10 




ESO488049 



P[M/pc 3 | 

(l ill _ 




FIG. 9: Density profile fits for the galaxies ESO1870510, 
ESO0488049, ESO0840411. 
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TABLE I: Galactic Parameter Values 
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RESUMEN 

El campo escalar con un potencial dc cosh sc comporta exactamente del mismo 
modo que la materia obscura fna (CDM), en la region donde el campo escalar 
oscila alrcdcdor de su mfnimo. Del mismo modo, en el regimen lineal, la hipotcsis 
de materia obscura como campo escalar (SFDM) predice la misma formacion de 
estructura que la predicha por la materia obscura fna. Estos rcsultados implican 
que los modclos CDM y SFDM son cquivalentcs desde el punto de vista cosmologico. 
Los parametros libres del model SFDM quedan fijos por medio de las observaciones 
cosmologicas. En un trabajo anterior mostramos que al resolver las ecuaciones de 
Einstein Klein Gordon utilizando dichos valores para los parametros libres, el campo 
escalar se colapsa y forma objetos estables con masa del orden de 10 12 M Q . En el 
presente trabajo utilizamos soluciones anahticas para las ecuaciones de Einstein 
Klein Gordon en el caso piano y en el Hmite de campo debil. Con estas soluciones, 
mostramos que el pcrfil de densidad del campo escalar corresponde al de un halo 
con una densidad central casi plana y que dicho halo coincide con cl del modclo 
de CDM en una amplia zona de la region exterior. Dicho resultado podrfa resolver 
el problema de los picos dc densidad en los halos de materia obscura, sin ninguna 
hipotesis adicional, lo que apoya la viabilidad del modelo de SFDM. De ese modo, 
el modelo de SFDM se puede considerar como un modelo alternative al de CDM 
con un WIMP auto-interactuante. 

ABSTRACT 

In the region where the scalar field endowed with a cosh potential oscillates 
around its minimum, it behaves exactly in the same way as the cold dark matter 
(CDM) model. Also, in the linear regime, the scalar field dark matter (SFDM) 
hypothesis predicts the same structure formation as the cold dark matter model. 
This means that CDM and SFDM are equivalent from the cosmological point of 
view. The free parameters of the SFDM model can be determined by means of cos- 
mological observations. In previous work we showed, by solving the Einstein Klein 
Gordon equations, that if we use such parameters, the scalar field collapses forming 
stable objects with a mass around 10 12 Af Q . In the present work we use analytical 
solutions of the flat and weak field limit of the Einstein- Klein-Gordon equations. 
With these solutions we show that the scalar field density profile corresponds to a 
halo with an almost flat central density and that this halo coincides with the CDM 
model in a broad outer region. Such a result could solve the problem of the density 
cusp DM halo in galaxies without any additional hypothesis, thus supporting the 
viability of the SFDM model. In this way, the SFDM model can be consider as an 
alternative model to the CDM one, with a self interacting WIMP. 

Key Words: COSMOLOGY: DARK MATTER GALAXIES: HALOS 
GALAXIES: DENSITY PROFILES OSCILLATON 
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1. INTRODUCTION 

The Lambda Cold Dark Matter (ACDM) model 
has recently shown an enormous predictive power. 
It can explain the structure formation of the Uni- 
verse, its accelerated expansion, the micro Kelvin 
fluctuation of the Cosmic Microwave Background 
Radiation, etc. Nevertheless, some issues around 
this model related to the formation of galaxies have 
arisen since the time it was originally proposed and 
remain to date. The CDM paradigm predicts a den- 
sity profile which corresponds to the Navarro-Frenk- 
White (NFW) profile, Navarro et al, 1997, given by 

PNFW = ji/ji T i\2 " W 

However, this profile seems to have some differences 
to the observed profiles of LSB galaxies. In this work 
we show that a flat central profile naturally arises 
within the scalar field dark matter hypothesis, im- 
plying that the central region of galaxies can distin- 
guish between CDM and SFDM. 

We work within the specific context of the so- 
called 'strong, self-interacting scalar field dark mat- 
ter' (SFDM) hypothesis that has been developed by 
several authors, Guzman and Matos 2000, Matos and 
Urena-Lopez 2000, 2001; Urena-Lopez, Matos and 
Becerril 2002; Matos and Guzman 2001; Alcubierrc 
ct al. 2002, 2003; Urena-Lopez 2002; Bohmer and 
Harko 2007; (see also Peebles 2000). A first pro- 
posal of the SFDM hypothesis appeared in a couple 
of papers by Ji and Sin 1994. They took a massive 
scalar field and were able to fit observations com- 
ming from some galaxies, taken into account also the 
contribution of baryons. A key point in this work 
was the use of the so-called excited configurations 
in which the radial profile of the scalar fields has 
nodes. From that, they determined that the mass of 
the scalar field should be of order of ~ 10~ 24 cV. 

A next proposal appeared in a paper by Schunck 
1998. He shows that a massless complex scalar field 
can be used as a dark matter model in galaxies to 
fit the rotation curves. In this model the internal 
frequency of the field plays the role of an adjustable 
parameter, and the radial profile of the scalar field 
also has nodes. 

However, as pointed out in Guzman and Urena- 
Lopez 2005 (see also Guzman and Urena-Lopez 
2004), this last proposal cannot be realistic because 
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a massless scalar field (whether real or complex as in 
Schunck 1998) cannot form a gravitationally bound 
configuration (see Seidel and Suen 1994). 

The key idea of the SFDM scenario is that the 
dark matter responsible for structure formation in 
the Universe is a real scalar field, <f>, minimally cou- 
pled to Einstein gravity with self-interaction para- 
metrized by a potential energy of the form (see also 
Sahni and Wang 2000) 

F($) = V a [cosh(A^$) - 1] , (2) 

where Vq and A are the only two free parameters of 
the model, kq — 8irG and we employ natural units 
Ti = c = 1. The effective mass of the scalar field is 
given by m| = KoV^A 2 . 

The advantage of the SFDM model is that it is 
insensitive to initial conditions and the scalar field 
behaves as CDM once it begins to oscillate around 
the minimum of its potential. In this case, it can 
be shown (see Matos and Urena-Lopez 2000, 2001) 
that the SFDM model is able to reproduce all the 
successes of the standard ACDM model above galac- 
tic scales. 

Furthermore, it predicts a sharp cut-off in the 
mass power spectrum due to its quadratic nature, 
thus explaining the observed dearth of dwarf galax- 
ies, in contrast with the possible excess predicted 
by high resolution N-body simulations with standard 
CDM, see Matos and Urena-Lopez 2001. 

The best-fit model to the cosmological data can 
be deduced from the current densities of dark matter 
and radiation in the Universe and from the cut-off in 
the mass power spectrum that constrains the number 
of dwarf galaxies in clusters. The favored values for 
the two free parameters of the scalar field potential 
are found to be, Matos and Urena-Lopez 2001: 

A ~ 20, 

V ~ (3x lQ- 27 m Pl )\ (3) 

where mpi = G -1 / 2 f=a 10 _5 g is the Planck mass. 

This implies that the effective mass of the scalar 
field should be to* ~ 9.1 x 10~ 52 m Pl = 1.1 x 10~ 23 
eV. 

Let us explain why we suspect that the scalar 
field could be the dark matter at galactic scales as 
well. There are three main reasons. 

The first reason is that numerical simulations 
suggest that the critical mass for the case considered 
here, using the scalar potential (2), and the parame- 
ters given by Eq. (3), is approximately, Alcubierre 



et. al. 2002 



M r , 



O .l^L=2.5xlO 13 M . (4) 



This was a surprising result. The critical mass of 
the model shown in Matos and Ureha-Lopez 2000, 
2001, is of the same order of magnitude of the dark 
matter content of a standard galactic halo. Observe 
that the parameters of the model, Eq. (3), were fixed 
using cosmological observations. The surprising re- 
sult consisted in the fact that using the same scalar 
field for explaining the dark matter at cosmological 
scales, it will always collapse with a preferred mass 
which corresponds to the halo of a real galaxy. Thus, 
this result is a prediction of the cosmological SFDM 
model for galaxy formation. 

The second reason is that during the linear 
regime of cosmological fluctuations, the scalar field 
and a dust fluid, like CDM, behave in the same way. 
The density contrast in CDM and in the SFDM mod- 
els evolve in exactly the same form and then both 
models predict the same large scale structure for- 
mation in the Universe (see Matos and Urena-Lopcz 
2001). The differences between the CDM and SFDM 
models begin to appear in the non linear regime of 
structure formation, so that there will be differences 
in their predictions on galaxy formation. 

The third reason is the topic of this work. A 
scalar field object (e.g. an oscillaton) contains a flat 
central density profile, as seems to be the case in 
galaxies. 

In the case of the SFDM, the strong self- 
interaction of the scalar field results in the formation 
of solitonic objects called 'oscillatons', which have a 
mass of the order of a galaxy (see for example Urefia- 
Lopez 2002, Urena-Lopez et al 2002 and Alcubicrre 
et al 2003. Also Seidel and Suen 1991, 1994, Hawley 
and Choptuik 2000, and Honda and Choptuik 2001). 
In this work we will show that these models contain 
an almost flat central density profile, i.e., they do 
not exhibit the cusp density profiles characteristic of 
the standard CDM hypothesis. 

Before starting with the description, we want to 
emphazise the fact that the scalar field has no in- 
teraction with the rest of the matter, thus, it does 
not follow the standard lines of reasoning for the 
particlc-likc candidates for dark matter. The scalar 
field was not thermalized, that is, the scalar field 
forms a Bose condensate, and thus behaves strictly 
as cold dark matter from the beginning. 

The rest of the paper is organized as follows. In 
the next section we use the fact that Galaxies have 
a weak gravitational field and thus their space-time 



is almost flat. The main goal of this section is to 
study the physics provoking the flatness behavior 
of the density profiles at the center of the oscilla- 
tons. Some results of this section intersect with those 
presented by Jin and Sin 1994, where they studied 
the behavior of the weak field limit of a complex 
scalar field. We remark that we do our analysis for 
a real scalar field and some differences do arise due 
to a different current conservation. In section 3 we 
present our contribution making an analysis of the 
Einstein-Klein-Gordon (EKG) equations in the rela- 
tivistic weak field limit and solve them for the per- 
turbed metric coefficients. Then, we compare these 
solutions with the ones obtained by solving numeri- 
cally the complete EKG system and using the whole 
potential (2), we show that the relativistic weak field 
limit is indeed a very good approximation. In section 
4 we compute the energy density of the scalar field 
obtained in section 3, and compare it with actual 
observations of LSB galaxies from which the den- 
sity is inferred from the rotational curves, showing a 
good matching in the external regions and a match 
at least similar in some of the internal regions, in 
any case, better than a fit with a cusp-like behav- 
ior of the density. Finally in section 5 we give our 
conclusions. 

2. PHYSICS OF THE SCALAR FIELD. FLAT 
SPACE-TIME CASE 

In this section we derive the physics of the scalar 
field in an analogous way as it was done for the com- 
plex scalar field by Jin and Sin 1994 (see also Lee and 
Koh 1996). In a normal dust collapse, as for example 
in CDM, there is in principle nothing to avoid that 
the dust matter collapses all the time. There is only a 
radial gravitational force that provokes the collapse, 
and to stop it, one needs to invoke some virializa- 
tion phenomenon. In the scalar field paradigm the 
collapse is different. The energy momentum tensor 
of the scalar field is 



2V(*)] . (5) 



We will consider spherical symmetry, and work with 
the line clement 



ds 2 



J dt 2 + e 2 ^dr 2 + r 2 dif , 



(6) 



with \i = /j,(r, t) and v = v (r, t), being this last func- 
tion the Newtonian potential. The energy momen- 
tum tensor of the scalar field has then the compo- 
nents 



— T°o — p<t> — 



1 



'$ 2 + e" 2 ^' 2 + 2V($\ 7) 



3 



Toi = r<s, = 

T 1 = Pr = K 



T 2 2 = PX = 



1 



(8) 

e -2 M$ /2 _ 2^($)1 (9) 
-2^<j>2 _ g-2^^/2 _ 2V{<§>) 



(10) 



and also T 3 3 = T 2 2- These different components are 
identified as the energy density /?$, the momentum 
density V$, the radial pressure p r and the angular 
pressure p±. The integrated mass is defined by 



M{x) 



Air [ p<t,(X)X 2 dX . 
Jo 



(11) 



The radial and angular pressures are two nat- 
ural components of the scalar field which stop the 
collapse, avoiding the cusp density profiles in the 
centers of the collapsed objects. This is the main 
difference between the normal dust collapse and the 
SFDM one. The pressures play an important roll 
in the SFDM equilibrium. In order to see this, and 
considering that galaxies are almost flat, we conclude 
that the Newtonian approximation should be suffi- 
cient to describe the processes. In this section we 
will take the flat space-time approximation. 

Thus, we study a massive oscillaton without 
self-interaction (i.e. with potential V = lm|$ 2 ), 
in the Minkowski background (/i ~ v ~ 0). Even 
though it is not a solution to the Einstein equations 
as we are neglecting the gravitational force provoked 
by the scalar field, the solution is analytic and it 
helps us to understand some features that appear in 
the non-flat oscillatons. 

In a spherically symmetric space-time, the Klein- 
Gordon equation if 1 ® d a dfi<& — dV/d& = 0, (where 
r] a "d a dj3 stands for the D'Alambertian), reads 



(12) 



where over-dot denotes d/dt and prime denotes 
d/dr. The exact general solution for the scalar field 
$ is 



*(i,r) 



D ±ikr 



(13) 



and we obtain the dispersion relation fc 2 = u> 2 — m\ . 
For lu > m$ the solution is non-singular and vanishes 
at infinity. We will restrict ourselves to this case. 
It is more convenient to use trigonometric functions 
and to write the particular solution in the form 

$(t,a;) = $o^^cos(art), (14) 



where x = kr. It oscillates in harmonic manner in 
time. The scalar field can be considered to be con- 
fined to a finite region, see Urena-Lopez 2002, and 
Urena-Lopez et al 2002. 
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Fig. 1. The energy density profile for the Scalar Field 
Dark Matter model (see also Schunck 1997). The para- 
meters used in this plot are m$ = 20, u> = 20.001, $o = 
5 10" 2 . The density is given in arbitrary units and r is 
given in kpc. 

The analytic expression for the scalar field energy 
density derived from Eq. (14) is 



Vfc 2 

2x 2 



sin(a;) 



cos(x) 



(15) 



\x)) cos 2 (wt) +uj 2 sin 2 (x)] 



which oscillates with a frequency 2u>t. Observe that 
close to the central regions of the object, the density 
of the oscillaton behaves like 



1 



$ 2 /fc 2 [u 



k 2 cos 2 



M)]+0(.t 2 ), (16) 



which implies that when x — * the central density 
oscillates around a fixed value. 

On the other hand, the asymptotic behavior 
when x — * oo, is such that p$ ~ l/x 2 , i.e. far away 
from the center, in this approximation, the flat oscil- 
laton density profile behaves like the isothermal one. 
The mass function oscillates around M ~ usual 
for the galactic halos. 

In order to understand what is happening within 
the object, observe that the KG equation can be 
rewritten in a more convenient form in terms of the 
energy density, as 



dp<t> 
dt 



1_ d_ 
r 2 dr 



2 V<j,) = 0. 



(17) 



This last equation has a clear interpretation: 
Since its form looks like the conservation equation, 



p + V • J = 0, equation (17) represents the conser- 
vation of the scalar field energy. It also tells us that 
there is a scalar field current given by 



J$ = —T<s>r 



sin(x) sin(2cjt) _ 
a;cos(a;j — sm(a;)J 5 r. 



Observe that the quantity involved in this current 
is the scalar field momentum density (8). Although 
the flux of scalar radiation at large distances does 
not vanish, there is not a net flux of energy, as it can 
be seen by averaging the scalar current on a period of 
a scalar oscillation. We also see that the only trans- 
formation process is that of the scalar field energy 
density into the momentum density, and viceversa. 
For the realistic values (3) this transfer is very small. 

In Fig. 1 we show the behavior of the SFDM den- 
sity profile for a typical galaxy and in Fig. 2 we show 
the comparison between the NFW, the isothermal 
and the SFDM density profiles for the same galaxy. 
Observe that the SFDM and NFW profiles remain 
very similar up to 10 kpc, then the SFDM profile 
starts to follow the isothermal one. 
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Fig. 2. Comparison between the energy density profile for 
the Scalar Field Dark Matter model with the NFW and 
the Isothermal models. The parameters for the isother- 
mal model are p Iso = 0.3/(r 2 + 8 2 ) and for the NFW 
profile are Pnfw = W/(r(r + 8) 2 ). The parameter used 
for the SFDM model are the same as in the previous 
figure. 

In the next section we will see that if the 
gravitational force is taken into account, the oscil- 
laton is more confined (see also Urena-Lopcz 2002, 
Ureha-Lopez et al 2002, and Alcubierre et al 2003). 
The parameters used in the figures, correspond to a 
middle size galaxy. 



3. WEAK FIELD LIMIT EQUATIONS 

In this section we derive a novel method for in- 
tegrating the perturbed EKG equations and show 
that the solutions are in very good agreement with 
the numerical ones. This allow us to use these so- 
lutions to fit the rotation curves of several observed 
LSB galaxies. 

Within general relativity, the evolution of the 
scalar and gravitational fields are governed by the 
coupled EKG equations, the last one appearing from 
the conservation of the energy-momentum tensor 



rpCtfi 



KoSaf3, (18) 

$< q (V 2 -to 2 )$ = 0, (19) 



here R a p is the Ricci tensor, S a p = T a p — 
(l/2)g a0 T\, and V 2 = {VV^g)d ll [y/^gg^d v \ is 
the covariant DAlambcrtian operator. 

For simplicity, we continue to consider the non- 
static spherically symmetric case, given by Eq. (6). 
As usual in the weak field limit, we suppose the met- 
ric to be close to the Minkowski metric 770.^3 



g a p = flap 



1af3, 



(20) 



where \h a p\ -C 1, then we will consider an expansion 
of the functions in the metric of the form 



e 2u{t,r) =1 + ( ?2V(t,r)+0(e i ), 
e 2„(t,r) = i + e ^ 2 U(t,r)+0(e 4 ), 



(21) 



where e is an expansion parameter. We also consider 
that the spatial and time derivatives of the geometric 
quantities are regarded like 



d_ d_ 
dt dr 



(22) 



then to first order in e 2 the Ricci tensor components 
are respectively 



[V :rr -U, tt H — V, r ], 
r 

2 

[U, tt — V, rr + — U, r ], 



Rtr = [~U, t ]. (23) 
r 

On the other hand, the source is computed in the 
flat space in this case as well, thus the scalar field 
satisfies to e 2 order the Eq. (12) as in the previous 
section. 

It is important to emphasize that the relation 
(22) is the lowest one in the geometric fields, but 
it does not consider small velocities for the sources. 
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This is different from the Newtonian limit where the 
derivative relation for the scalar field is d r ~ edt and 
d t ~ ed r for the geometric fields (see Seidel & Suen 
1990, and Guzman & Ureha-Lopcz 2004). 

Consistent with the T a p computed in the flat 
space-time, the right hand side elements in Einstein's 
equations, are written as 



S tt = $,t $,t 



-m 1$ 



S tr = $, r $, t . 



1 

1 



2 I -=is 1 2 



$, r $, r +-m 2 |$| 2 , 



(24) 



In this case it is convenient to introduce the dimen- 
sionlcss quantities 



t = mt, Q, = 



ui 
m 



(25) 



where we note that the bosonic mass m is the natural 
scale for time and distance. In terms of these new 
variables the general solution to Eq. (13) takes the 
form 



$(r, x) — — cxp(±ixy/ Vt 2 — 1) cxp(±i£!t) 



(26) 



The physical properties of the solution depend on 
the ratio fl = uj/m. For VL < 1 the solution decays 
exponentially but it is singular at r = 0. On the 
other hand, f2 > 1 allows for non-singular solutions 
which vanish at infinity. We will restrict ourselves 
to this case. We will write the particular solution in 
the form 



(27) 



v / Ko$(r, x) = <j>(x) cos(Qt), 

where the spatial function is given by 

,. sin(xv / ^ 2 - 1) 
4>{x) = 0o • 



(28) 



Because of the functional form of the scalar field 
(27) we introduce the following ansatz for the metric 
perturbations (21) 

V(t,x) = V (x) + V 2 (x)cos(2nT), 
U(t,x) = U (x) + U 2 (x)cos{2Qt). (29) 

We adopt the following method to solve the EKG 
equations. In the first approximation we substitute 
the solution of a lower approximation (the flat case) 
into the next approximation (with e 2 ) and solve the 
resulting differential equations. As we will see, this 
standard approximation works well in our case. In 



terms of these expressions the Einstein's equations 
Ra0 = KoS a p finally read 



-U 2 



x 



V 2 , xx +-V 2 , x +4ft 2 £/ 2 
x 



-Vq, X x + — Uq, x 
X 



-V 2 , XX +-U 2 , X -AQ 2 U 2 
x 



2 V 2 
2 V 2 



(30) 



3.1. Scaling properties 

^From system (30) we know that the scalar field's 
maximum amplitude 0(0) = (po^/fl 2 — 1 could be 
taken as the expansion parameter e and in this case 
£1 must be of order 1. Then it is always possible to 
solve the system (30) ignoring e and replacing <j) by 
its normalized function 



4>(x) 



sin(xVtt 2 - 1) 
xy/Vt 2 - 1 



(31) 



Solutions <j), Uq, U 2 , Vq, V 2 of this normalized system 
depend only of the arbitrary characteristic frequency 
which modulate the wave length of </>. On the 
other hand, for each value of f2 there is a complete 
family of solutions of the scalar field <fi and the metric 
perturbations h a p which are related to each other by 
a scaling transformation characterized by 4>o 



htt = 



^oVfi 2 - l<£cos(fir), (32) 
(j)l{tt 2 - 1)[2J7 + 2J7 2 cos(20r)], 
-cfl{n 2 - \)[2V Q + 2V 2 cos(2Qr)], 



In this context the weak field limit condition h a /3 *C 
1 translates into 

^l(fl 2 -l)\2V\ < 1, & cj) 2 (n 2 ~l)\2U\ < 1. (33) 

Here we will introduce a specific notation for the 
spatial functions of the metric perturbations: 





4(n 2 


- 1)2U , 


(34) 






- 1)2U 2 , 






0o 


- 1)2V , 




h {2) - 


2 (ft 2 


- l)2V 2 . 





3.2. Metric perturbations solutions 

The system of equations Eqs. (30) can be solved 
and the spatial functions of the metric perturbations 
have analytic solutions given by 



U 2 = 



1 



8(fi 2 - 1) 



sm 2 (xV^ 2 - 1) 



+ 



Vn 2 - 1 sm(2x^/n 2 - 1) 



v tt 



2 x 
2 1) rsin(2a;v / ^ 2 - 1) 



(20 
8(0 2 - 1 



2xV^ 2 - 1 



Ci{2xy/W - 1) + 



+ ln(2x\/0 2 - 1; 



C 



voi 



C 



V02, 



Vo 



8(^2 



V^ 2 - 1 sin(2xV0 2 - 1) 
2 x 



+ Ci(2aVn 2 - 1) - 

- ln(2aVfi 2 - 1)1 + Cy 22 , 



1 



8(fi 2 - l 

1 sin(2xV0 2 - 1) 
~ 2 



1 1 cos(2x v / r2 2 - 1) 
2x2 + 2~ 



xV^ 2 - 1 



V01 



+ CW,(35) 



where C« is the cosine integral function and Cyoi> 
CV02, CV22, and CVoi are integration constants. 

3.3. Weak Field Validity Range 

^From equation (13) it is evident that, in the 
limit in which we are working with, the KG equation 
is decoupled from Einstein equations. Imposing reg- 
ularity at the origin and asymptotic flatness to the 
KG solution we have chosen (27) with (28) as our 
scalar field particular solution where <p and fl > 1 
are still free parameters. 

On the other hand, regularity at x = requires 
h rr (x = 0,t) = which implies 



CVoi = 0, Cxjqx = 



o 2 



8(n 2 - i)' 



(36) 



then for the perturbations Cyo2 and CV22 are still 
free integration constants. Now we will describe the 
asymptotic behavior of these perturbations. Due to 
U2 being at least one order of magnitude smaller 
than Uq and its value oscillating around zero, it is Uq 
which determines the behavior of h rr . The Uq value 
starts to oscillate, very near to the origin, around 
Cf/oi keeping this behavior asymptotically. Then the 
asymptotic value of h rr is the finite Cjjoi value. Con- 
trary to this h u , due to the logarithm terms in Vo 



and V 2 , htt, is singular at infinity. Thus, the weak 
field condition (33) is fulfilled only in a finite spatial 
region around the origin, i.e., due the the approxi- 
mation the solution is contained in a box, for which 
the walls are sufficiently far away from the center 
of the solution. We will say that this is the region 
where our weak field approximation is valid. 

Unique solutions for the EKG system will be ob- 
tained fixing the <pQ and parameters and the con- 
stants CV02 and CV01 within the validity range of the 
approximation. As it is known the potentials mea- 
surement does not have physical sense by themselves, 
then unique solutions will be determined through 
metric dependent observable quantities. Using the 
expressions given by Eqs. (34), we can obtain the 
perturbed metric functions in terms of these solu- 
tions. In Fig. 3 we present a plot of these metric 
perturbation functions, as well as of the scalar field, 
for two values of £1. 

n=i.4 




£1=1.1 




Fig. 3. Profiles of <j>(x), h„(x), h„(x),h^(x) and 
hip (x) with cj>o = 1; see text for details. 
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It is important to mention that the width of the 
validity region, where (33) is fulfilled, depends on SI 
and 4>o. This is because it is the factor 0g which mod- 
ulate the perturbations, (see (32)). What it is not 
evident until the solutions (35) are observed is that 
the SI value, independently from <f>Q, could make the 
validity range width bigger. This is because in the 
logarithm argument there is the expression y/fl 2 — 1, 
then as SI is closer to one the logarithm terms rise 
more slowly. 

The order of magnitude for the other parameter 
(fro in the metric perturbations, can be naturally de- 
termined from the asymptotic value taken by h rri 
which is reached very close to the origin 



lim \hr 



4 V 



(37) 



As SI is nearly 1 the h rr magnitude is given by cV 
It is well known that for week field systems like our 
Solar System the metric perturbations go like h a ~ 
10 -6 . This value restricts our d>n to be <bn < 10 -3 . 
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Fig. 4. (Left) Spatial function cj> of the scalar field 
$. Solid and da sh lines are analytical solutions 
with (f> = 0.001/^(S1 2 - 1) = 0.001591, and <j> = 
0.002/^(^2 - 1) = 0.001946, respectively. The crosses 
are the corresponding numerical solutions with <^>i(0) = 
0.001, n = 1.181008 and <£i(0) = 0.002, fi = 1.433822; 
see text for details. (Right) Spatial functions /ijt and 
h u of the metric perturbations. 



3.4. Analytical Solutions vs Numerical Solutions 

Analytic solutions <f>o, h^J(x), h$)(x), (x) 
and hi® (x) are shown in Fig. 3. The value of tp 
is 1 in both plots. As already was noticed, the value 
of SI characterizes each family of solutions. Mainly 
SI determines the wave length of <f> and the increase 
rate of and ; as SI is closer to 1 this rate is 
smaller. These characteristics are shown in Fig. 3. 

The exact EKG equations in spherical symme- 
try with a quadratic potential, were solved numer- 
ically in Ureha-Lopez 2002, and Urena-Lopez et al 
2002) and found the so called oscillatons. In those 
works boundary conditions are determined by requir- 
ing non-singular and asymptotically flat solutions, 
for which the EKG become an eigenvalue problem. 
The free eigenvalue is the scalar field's central value 
4>i(x = 0) which labels the particular equilibrium 
configuration, and the fundamental frequency SI is 
an output value. In those works it is also noted that 
weak gravity fields are produced by oscillatons with 
<pi(x = 0) <C 1. In Fig. 4 we compare some of these 
numerical solutions (NS) with the analytical solu- 
tions (AS) within a central region. The constant 
values of the AS are fixed to fit better the NS in- 
side the weak field validity range. From these plots 
we can conclude that our solutions are a very good 
approximation for the exact EKG equations in the 
weak field limit. The principal advantage of this ap- 
proximation is the analytical description of the solu- 
tions. 




X 



Fig. 5. (Left) Energy-Momentum density profiles for two 
scalar field configurations one with (f>o = 1.159f f 0~ 3 and 
n = 1.181008 and the other with <j> = 5.2225 10" 3 and 
SI = 1.153246. (Right) Difference between density pro- 
files from the exact EKG equations (pjv) and the density 
computed in the flat space-time (p); see text for details. 



4. SCALAR FIELD AS DARK MATTER: HALO 
DENSITY PROFILE 

In this section we explore whether or not the 
scalar field could account as the galactic DM halos. 
Specifically wc compare the SFDM model density 
profile and the profiles inferred throughout the ro- 
tation curves profiles of galaxies which are mostly 
formed by DM. 

As long as we are concerned with perturbations of 
the flat space-time due to the scalar field, we do not 
consider the baryonic matter gravitational effects, we 
expect that our approximation will be better suited 
for galaxies with very small baryonic component. 

We will compare the energy-momentum density 
for the scalar field given by Eq. (7), in the relativis- 
tic weak limit approximation, where for the metric 
functions, Eq.(21), we use the solutions to the per- 
turbations given by Eqs. (29,35). 



ES OO 1 40040 ES OO 1 40040 




Fig. 6. Density profile fits for the galaxies ESO0140040, 
ESO2060140, the right panel is a zoom of its central 
region, it is plotted in order to visualize complete er- 
ror bars. The horizontal line is in kpc, vertical line in 
Mq/pc . See text and Table 1 for fit details. 

This is consistent with the fact that gravity does 
not modify the scalar field behavior. This approxi- 
mation in the weak gravitational field limit is very 
good as we can see in Fig. 5. In those plots we 
show the energy-momentum density profiles for two 
scalar field configurations with different maximum 
amplitudes at the origin (f)(0) — O \/^ 2 — 1- It is 
important to notice that as (f)(0) decreases the gravi- 
tational field gets weaker, and the difference between 
the density from the complete EKG equations and 
from our approximation is smaller. 

The density profile fits allow us to obtain an es- 
timation of the parameters at the galactic level: the 
fundamental frequency O and the scalar field con- 
stant 4>q. The third parameter involved in the den- 
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Fig. 7. The same as Fig. 6, but for the galaxies 
ESO3020120, ESO3050090. 

sity profiles is the scalar field mass, we will fix it to be 
m = 1CT 23 eV. This value was fitted for the SFDM 
model from cosmological observations in Matos & 
Urena-Lopez 200 f. 




Fig. 8. The same as Fig. 6, but for the galaxies 
ESO4250180, ESO1200211. 

4.1. Density Profile fits 

The first qualitative feature of the energy-density 
profile that we want to emphasize is that in the cen- 
tral region it is non cuspy (see Fig. 5). It is im- 
portant to take into account that instead of density 
profiles, rotation curves are the direct observable for 
galaxies. Nevertheless, for galaxies dominated by 
DM, their rotation curves could model the DM den- 
sity profile more trustfully. We choose a subset of 
galaxies from the set presented in McGaugh 2001, 
the common characteristic for the selected galaxies 
is that the luminous matter velocity contribution to 
the rotation curves is almost null. 



With the scalar field mass m fixed, the profiles 
fits were made for the Q and <f>o values with good 
X 2 statistic, see Table 1. In most of the cases the 
non central observational data were the good fitted 
points, those data points also have smaller error bars. 
The density profile fits are in Figures 6,7,8,9, where 
we show several galaxies with the density computed 
from the observed rotational profile versus the den- 
sity obtained with our SFDM description. In some 
of them we were able to compare regions within less 
than 0.5 Kpc. 

In Table 1 the fundamental frequency fi is listed 
for each galaxy. We found that the temporal de- 
pendence for the energy-momentum density profile 
is harmonic with a temporal period T = n/Cl. The 
column Ap(0) corresponds to the maximum change 
in the central density for a period of time T . Finally 
for all the galaxies the 4>o value is well inside of the 
weak gravitation field limit <fio < 10~ 3 . 

5. CONCLUSIONS AND FUTURE PROSPECTS 

We have found analytic solutions for the EKG 
equations, for the case when the scalar field is con- 
sider as a test field in a Minkowski background, and 
in the relativistic weak gravitational field limit at 
first order in the metric perturbations. With these 
solutions we have shown that non-trivial local behav- 
ior of the scalar field holds the collapse of an object 
formed from scalar field matter. The scalar field con- 
tains non trivial, natural effective pressures which 
stop the collapse and prevent the centers of these 
objects from having cusp-Ike density profiles. Even 
within this simple approximation it has been possi- 
ble to fit, with relative success, the density profiles 
for some galaxies showing non cuspy profiles. 

Together, all the features of the SFDM model al- 
low one to consider this model as a robust alternative 
candidate to be the dark matter of the Universe, as 
was suggested by Guzman and Matos 2000, Matos 
and Guzman 2000, 2001, and Matos et al 2000. Fur- 
thermore, it has been shown previously that dark ha- 
los of galaxies could be scalar solitonic objects, even 
in the presence of baryonic matter, Hu et al 2000; Lee 
and Koh 1996; Arbey et al 2001, 2002; Sin 1994; and 
Ji and Sin 1994. Actually, the boson mass estimated 
in all these different approaches roughly coincides 
with the value to$ ~ 10 _23 eU, even if the later was 
estimated from a cosmological point of view, Matos 
and Ureha-Lopez 2001. We can appreciate the non- 
trivial characteristics of the proposed potential (2): 
Its strong self-interaction provides a reliable cosmo- 
logical scenario, while at the same time it has the 
desired properties of a quadratic potential. Finally, 



TABLE 1 
GALACTIC PARAMETER VALUES 



Galaxy 




Q 




00 




x 2 


p(r = 0) 


Ap(0) 


T 
















[M Q /pcs 3 ] 


[M Q /pc S 3 } 


[yrs] 


ESO0140040 


1 - 


f-8 10" 


-9 


1.87 10" 


-3 


12.366 


0.569 10- 1 


0.755 10~ 8 


1.603382750 10 7 


ESO0840411 


1 - 


f 6 10- 


9 


5.95 10" 


-4 


1.338 


0.433 10~ 2 


0.147 10~ 9 


1.603953416 10 7 


ESO1200211 


1 4 


-53 10 


-9 


2.04 10" 


-4 


10.062 


0.448 10~ 2 


0.530 10~ 9 


1.603382679 10 7 


ESO1870510 


1 4 


- 12 10 


-9 


3.28 10" 


-4 


3.190 


0.265 10- 1 


0.699 10~ 8 


1.603382570 10 7 


ESO2060140 


1 4 


- 18 10 


-9 


9.18 10" 


-4 


65.421 


0.308 10- 1 


0.206 10~ 8 


1.603382735 10 7 


ESO3020120 


1 - 


-29 10 


-9 


6.88 10- 


-4 


16.099 


0.279 10- 1 


0.170 10~ 8 


1.603382718 10 7 


ESO3050090 


1 4 


-40 10 


-9 


4.72 10- 


-4 


1.224 


0.181 10" 1 


0.243 10" 8 


1.603382699 10 7 


ESO4250180 


1 - 


f 4 10" 


9 


1.27 10" 


-3 


5.221 


0.132 KT 1 


0.105 10" 8 


1.603382757 10 7 


ESO4880049 


1 


f 3 10- 


9 


7.86 10" 


-4 


11.410 


0.377 10- 1 


0.212 10~ 8 


1.603382715 10 7 



the results presented here fill the gap between the 
successes at cosmological and galactic levels. 
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